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We present first-principles calculations of the coupling of quasiparticles to spin fluctuations in 
iron selenide and discuss which types of superconducting instabilities this coupling gives rise to. We 
find that strong antiferromagnetic stripe-phase spin fluctuations lead to large coupling constants 
for superconducting gaps with -symmetry, but these coupling constants are significantly reduced 
by other spin fluctuations with small wave vectors. An accurate description of this competition 
and an inclusion of band structure and Stoner parameter renormalization effects lead to a value 
of the coupling constant for an -symmetric gap which can produce a superconducting transition 
temperature consistent with experimental measurements. 
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Introduction .—The discovery of superconductivity in 
iron-based compounds with transition temperatures 
higher than 50 Kelvin in 2008 [T] has generated consid¬ 
erable interest in recent years and led to intense research 
activity. These materials consist of iron atoms with ad¬ 
ditional pnictogen or chalcogen atoms located above and 
below the plane of the iron atoms. Additional intercalat¬ 
ing layers give rise to various families of iron-based com¬ 
pounds, such as the fill-, 122-, Ill-families [2]. Iron se¬ 
lenide (FeSe) belongs to the structurally simplest family, 
the 11-family, which does not contain a spacer layer be¬ 
tween the iron-selenium layer and exhibits superconduc¬ 
tivity with transitions temperatures up to 37 Kelvin un¬ 
der pressure Recently, there has also been consider¬ 

able interest in the properties of FeSe monolayers, where 
superconductivity with transition temperatures exceed¬ 
ing 60 Kelvin has been reported 13 E). 

An important question in these materials is the nature 
of the microscopic pairing mechanism and the symme¬ 
try of the superconducting gap. As in the cuprates, su¬ 
perconductivity typically emerges in the iron-based com¬ 
pounds when an antiferromagnetic parent state is doped. 
This observation led to the proposal [7] that spin fluctu¬ 
ations (paramagnons) could act as the superconducting 
glue in these materials [7]. However, in contrast to the 
cuprates, the iron-based materials generally have multi¬ 
ple Fermi surfaces with several electron pockets at the 
M-point of the Brillouin zone (corresponding to the unit 
cell containing two iron atoms) and multiple hole pock¬ 
ets at the F-point. Mazin et al. suggested that spin- 
fluctuation mediated scattering of Cooper pairs between 
electron pockets and hole pockets gives rise to an s±- 
symmetric superconducting gap [7] , which has a constant 
absolute magnitude, but switches sign between electron 
and hole pockets. The sign-changing gap was predicted 
to give rise to a resonance in the neutron scattering spec¬ 
trum [8] , which was subsequently observed in several ex¬ 


periments laiin]. 

Many theoretical approaches have been developed to 
study spin-fluctuation mediated superconductivity in the 
iron-based compounds. In purely empirical approaches, 
both the electronic band structure and the interacting 
spin susceptibility are parametrized using experimental 
data (such as angle-resolved photoemission, nuclear mag¬ 
netic resonance and neutron scattering results) [TTl |T2]. 
In another approach usmi, a theoretical band struc¬ 
ture from a density-functional theory calculation is used 
to parametrize a tight-binding Hamiltonian with added 
interaction parameters (such as the Hubbard U or Hund’s 
J) adjusted to reproduce experimental findings. Then, 
superconducting properties of the resulting Hamiltonian 
are studied. 

While the aforementioned theories have been very in¬ 
structive, their applications have been limited by the 
availability of concrete experimental data needed to de¬ 
termine their input parameters. Hence, there is a need 
for a fully first-principles theory without empirical pa¬ 
rameters. We have recently developed such an ab ini¬ 
tio theory m for the spin fluctuation-electron coupling 
based on the work of Overhauser and coworkers on the 
homogeneous electron gas min]. 

In this paper, we apply our first-principles theory of 
the electron-spin fluctuation coupling to iron selenide, 
the structurally simplest iron-based superconductor. We 
reveal a complex interplay between different spin fluc¬ 
tuations. In agreement with experiment, superconduc¬ 
tivity depends sensitively on the height of the selenium 
atoms relative to the iron atom plane. In the vicinity of 
a transition to an antiferromagnetic stripe phase, super¬ 
conductivity with an -symmetric gap function becomes 
favorable. 

Methods .—The superconducting order parameter, the 
gap function, is typically expressed as A{k) = \A\g{k) 
with g{k) describing the symmetry of the gap function 
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[18] . To investigate which types of superconducting insta¬ 
bilities an effective spin-fluctuation-mediated electron- 
electron interaction Vsf gives rise to, we compute the 
electron-spin fluctuation coupling strength for different 
symmetry functions g{k) according to [18] 

^[9] = —DfX ( 1 ) 

Y^nk.n'k' 9 {k)VsF{nk, n'k')g{k')8{enk - ^F)8{en'k' - ^f) 
['Enk9‘^{k)S{enk-eFW 

( 2 ) 

where cf denotes the Fermi energy and Df is the den¬ 
sity of states per spin at the Fermi energy. We de¬ 
fine VsFi'^k^n'k') as the matrix element of the spin- 
fluctuation mediated interaction for the scattering of the 
spin-singlet Cooper pair {nk t, n—k |) to {n'k' t, n'—k' | 
) (with n denoting a band and k a /c-point in the Bril- 
louin zone). We follow the standard convention [18] and 
normalize g{k)^ such that Df = '^nk 9^(^)^(^nk — ^f)- 

The effective spin-fluctnation mediated interaction is 
obtained using a recently developed first-principles for¬ 
malism nann based on the work of Overhauser et al 
for the homogeneous electron gas [m HZ]. In this ap¬ 
proach, the effective interaction is expressed as the sum of 
a bare Coulomb interaction, a contribution arising from 
charge fluctuations and a contribution arising from spin 
fluctuations, which is given by 

VsF{r,r',uj) = 3x 

dri J dr2lxc{r,ri)xs{ri,r2,uj)I:cc{r2,r'), (3) 

with Ixc(r,r') = 5'^Exc/[dm(r)5m(r')\ {E^c is the 
exchange-correlation energy and m{r) the spin density). 
Also, Xs'(r, r', cj) denotes the interacting spin susceptibil¬ 
ity of a nonmagnetic system, which is obtained by solv¬ 
ing the Dyson-like equation of time-dependent density- 
functional theory m 

Xs{r,r',oj) = xo{r,r\uj) 

+ j dri J dr2Xo{r,ri,u})Ixc{ri,r2)xs{r2,r',u}), ( 4 ) 

with xo denoting the non-interacting susceptibility. 

The physically appealing effective spin fluctuation- 
mediated interaction of Overhauser et al. can be put 
on a firm theoretical footing by an analysis of Feynman 
diagrams mi HI- While this theory neglects certain dia¬ 
grams corresponding to nonlinear polarization processes, 
it should be valid in the vicinity of the Fermi surface and 
thus describe superconducting properties accurately. 

Computational details. —To obtain a mean-field theory 
starting point for the calculation of the spin-fluctuation 
mediated interaction, we carry out density-functional 
theory calculations in a plane-wave basis using the 
Quantum ESPRESSO program package [23]. We em¬ 
ploy the local-density approximation (Perdew-Zunger 


parametrization) and norm-conserving pseudopotentials 
with a 55 Ry energy cutoff. Eor iron, the non-linear core 
correction is used. We use experimental lattice constants 
(a = 3.77 A, c = 5.52 A) [^ of the tetragonal phase. 
We have carefully verified that our band structures and 
magnetic phase diagrams agree with all-electron results 
[251126]. 

We then compute the non-interacting susceptibilities 
on a 32 X 32 X 1 /c-point grid in the Brillouin zone us¬ 
ing 50 empty states and a plane-wave cutoff of 45 Ry. 
Eor this, we use the Berkeley GW program package m- 
Next, we compute the interacting spin susceptibilities. 
Within the local-density approximation, we need to eval¬ 
uate /^/^^(r,r') = f{n{r))6{r — r') with n(r) being 
the ground-state density (including the core contribu¬ 
tion) and f{n) is obtained analytically by computing the 
second derivative of the exchange-correlation energy with 
respect to the spin density. Einally, the Eermi surface 
averages in Eq. (§ are evaluated using a Lorentzian rep¬ 
resentation of the (5-function [Srj{uj) = I/tt x r]l{uj‘^ +^^)] 
with 77 = 0.1 eV. 

Results —. Eigure [^a) shows the magnetic moments 
of various magnetic phases of EeSe as a function of the 
selenium height, which has been identified as a crucial 
parameter for the occurence of superconductivity in this 
system and other iron-based compounds [28]. At low se¬ 
lenium heights, the system is nonmagnetic. As the height 
is increased, a transition to an antiferromagnetic stripe 
phase occurs. At even higher Se heights, an antiferromag¬ 
netic checkerboard phase also has a lower energy than the 
nonmagnetic state and a moment develops. Einally, at Se 
heights larger than 1.4 A a ferromagnetic moment is ob¬ 
served. Note that bulk EeSe is experimentally known to 
be nonmagnetic [3], while a recent ARPES experiment 
reported indications that multilayer EeSe is in a stripe 
phase state [29] . 

Eigure [^b) shows the DET-LDA band structure of 
EeSe at the experimental selenium height of 1.47 A [24]. 
We observe three hole pockets in the vicinity of the F- 
point and two electron pockets near the M-point (note 
that a Brillouin zone corresponding to a unit cell con¬ 
taining two iron atoms is used). While DFT-LDA band 
structures agree qualitatively with experimental ARPES 
measurements for many iron-based compounds, the ef¬ 
fective masses of the electron and hole pockets near the 
Eermi level are typically underestimated by a factor of 2 
or 3 [301 El]. A similar finding was reported in a recent 
ARPES experiment on EeSe multilayers [29]. In Eig.j^b) 
we also show the renormalized band structure obtained 
by dividing all DET-LDA energies by a factor of 3. We 
note that a recent Shubnikov-de Haas oscillation mea¬ 
surement reported significantly reduced Eermi surfaces 
[32] . Eurther experimental work is necessary to fully re¬ 
solve this issue. 

Eigure shows the real part of the DET-LDA static 
non-interacting and interacting spin susceptibilities at 
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FIG. 1. (a): Magnetic moment per iron atom in the stripe an¬ 
tiferromagnetic phase (red hlled dots), checkerboard antifer¬ 
romagnetic phase (blue empty dots) and ferromagnetic phase 
(magenta diamonds) as function of the height of the selenium 
atoms above (and below) the plane of the iron atoms, (b) 
DFT-LDA band structure of iron selenide at the experimen¬ 
tal selenium height (dashed blue lines) and the renormalized 
band structure (solid red lines), where all energies have been 
divided by a factor of 3 [29]. 


three Se heights in the vicinity of the transition from the 
nonmagnetic phase to the antiferromagnetic stripe phase. 
At the smallest Se height (1.10 A), we observe strong 
features in the non-interacting susceptibility near q = 0. 
In the interacting spin susceptibility, these features are 
strongly enhanced indicating that at these wave vectors 
the Stoner condition, which states that magnetism occurs 
if IxcXo{q) (in matrix notation) has eigenvalues equal to 
unity, is almost fulfilled. 

At a selenium height of 1.15 A [see Figs. |^) and c)], 
another feature atgr = (l/2, 1/2, 0)^ emerges in ad¬ 
dition to the structure near q = 0 in the non-interacting 
susceptibility. This feature results from spin fluctuations 
with stripe phase character. When interactions are in¬ 
cluded, this new feature, however, becomes weaker. 

Finally, at a selenium height of 1.20 A [see Figs. |^) 
and f)], the system is very close to the transition to the 
antiferromagnetic stripe phase, see Fig. [^a). Now, the 
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FIG. 2. Real parts of the static non-interacting [(a), (c), 
(e)] and interacting [(b), (d), (f)] spin susceptibilities of iron 
selenide for selenium heights of 1.10 A [(a) and (b)], 1.15 A [(c) 
and (d)] and 1.20 A [(e) and (f)]. Shown is the G = G' = 0 
component of the susceptiblity matrix in Fourier space. 


peak atgr = (l/2, 1/2, 0)^ becomes very strong and is 
further enhanced when interactions are included. Again, 
the divergence of the interacting spin susceptibility indi¬ 
cates that the Stoner criterion is almost fulfilled signaling 
the onset of a transition to a stripe phase. 

To gain further insight, we fit our results to a model 
where off-diagonal elements of the non-interacting sus¬ 
ceptibilities in a plane-wave basis are neglected (neglect¬ 
ing the so-called local-field effects). We then adjust the 
Stoner parameter Ixc, which now is a single number, un¬ 
til the model reproduces the diagonal elements of the 
previously computed interacting spin susceptibility. In 
this way, we extract Ixc = 0.43 eV (per Fe atom), which 
agrees very well with the Stoner parameter for bee iron 
Ixc = 0.46 eV [33|. In addition, we find that the same 
value of Ixc can be used for all selenium heights that were 
studied. This indicates that the Stoner parameter is only 
weakly dependent on the crystalline environment. 

Next, we compute the electron-spin fluctuation cou¬ 
pling strength and evaluate A[^] for the lowest order sym¬ 
metry functions g{k) of the tetragonal crystal. Specif¬ 
ically, we compute the coupling strength for s-wave 
[g{k) = 1], [g{k) oc cos{kx) ^ cos{ky)], 42_^2 [g{k) oc 

cos{kx) — cos{ky)] and dxy [g{k) oc sin{kx) sin(/c^)] sym¬ 
metries. Fig. [^a) shows our results for the various cou¬ 
pling strengths as function of the selenium height. In the 
singlet channel, the spin-fluctuation mediated interaction 
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FIG. 3. (a) Electron-spin fluctuation coupling constant for 
various gap symmetries as function of the selenium height. 
Positive coupling constants indicate that the system may un¬ 
dergo a transition to a superconducting phase at sufficiently 
low temperatures, (b) Contributions to the total electron- 
spin fluctuation coupling constant for an -symmetric gap 
from the gr-points in the Brillouin zone at a selenium height 
of 1.20 A. (c) Positive and negative contributions to X[s±] as 
function of the selenium height. 


between electrons is repulsive and therefore the coupling 
strength for s-wave symmetry is always negative. The 
vicinity to a magnetic phase transition gives rise to an al¬ 
most singular spin-fluctuation mediated interaction and 
large, negative coupling constants ^ —1. This re¬ 
pulsion arising from magnetic interactions is much larger 
than the repulsion from the screened Coulomb interac¬ 
tions in standard metals, where = —/i ~ —0.2 [34] . 


Note that this repulsion must also be overcome by other 
pairing mechanisms, such as phonons. 

Figure [^a) shows that d-wave coupling constants are 
quite small, but can be positive and thus give rise to su¬ 
perconductivity at sufficiently low temperatures. Most 
interestingly, we observe that the s± coupling constant is 
large and negative for small selenium heights, but quickly 
increases and eventually becomes positive near the tran¬ 
sition to a stripe phase reaching values of approximately 
unity. We expect that this value of the coupling con¬ 
stant is reduced by the other contributions to the effec¬ 
tive interaction (such as the bare Coulomb interaction 
and the charge fluctuation-mediated interaction) only by 
a small amount^ as these contributions are much more 
isotropic and thus produce small values when integrated 
against the anisotropic g{k) in Eq. (§• 

Figure, [^b) shows the contributions to the total s± 
coupling strength from all gr-points in the Brillouin zone 
at a selenium height of 1.20 A. We observe significant can¬ 
cellations between positive and negative contributions. 
Fig. [^c) shows the total positive and total negative con¬ 
tributions to the s± coupling strength as function of the 
selenium height. While the negative contributions re¬ 
main relatively constant, the positive contributions in¬ 
crease rapidly as stripe phase spin fluctuations become 
enhanced. These spin fluctuations scatter electrons from 
the hole pockets at the F-point to the electron pockets 
near the M-point. In the s± scenario, the superconduct¬ 
ing order paramater switches sign between the F- and 
the M-point and therefore such scattering events are fa¬ 
vorable for the emergence of superconductivity. On the 
other hand, spin fluctuations with wave vectors in the 
vicinity of the F-point scatter Cooper pairs only from 
electron pockets to other electron pockets or from hole 
pockets to other hole pockets. Their contribution to the 
s± coupling constant is negative. Suppression of such 
fluctuations (for example, by application of pressure) 
could provide a path towards higher transition tempera¬ 
tures mm- 

Finally, we discuss corrections to the presented first- 
principles framework that arise from the lack of self- 
consistency in our calculations. Specifically, the non¬ 
interacting spin susceptibility should be computed from 
a quasiparticle band structure including renormaliza¬ 
tion effects arising from many-electron interactions in¬ 
stead of the mean-field DFT-LDA band structure. While 
such renormalization effects are reproduced by “beyond- 
DFT” approaches, such as GW theory [35] or dynami¬ 
cal mean-field theory [25], good agreement with experi¬ 
mental ARPES measurements can be obtained by sim¬ 
ply dividing all DFT-LDA band energies by a factor a, 
i-e. ink = ^nk - CF ^ ink = ink/c^^ Based on a re¬ 
cent ARPES experiment on FeSe multilayers [29|, we use 
(a = 3, see Fig. Bb). This rescaling of the band energies 
renormalizes the non-interacting susceptibility according 
to Xo = ^Xo ^Iso the density of states at the Fermi 
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level D'p = aDp. 

In addition to the electronic band structure, also the 
Stoner parameter is renormalized in the vicinity of a mag¬ 
netic phase transition where spin fluctuations are sig¬ 
nificant [SilST]. To include this effect, we also rescale 
the Stoner parameter according to = IxdP- We use 
P = 4.17 which reproduces the experimental magnetic 
moment and critical doping strength in LaOFeAs [38] 
(note that ^ = a x s with s = 1.39 being the rescaling 
parameter for I^c neglecting band structure renormaliza¬ 
tion effects). 

We now compute the coupling constant for the s±- 
symmetric gap in FeSe at the experimental selenium 
height employing a renormalized band structure and 
Stoner parameter. This yields X[s±] = 0.30. Including 
the effect of charge fluctuations results in a slightly re¬ 
duced value of X[s±] =0.28. A crude estimate of the su¬ 
perconducting transition temperature Tc = 

(where we use ujsf = 15 meV, a typical energy scale 
for spin fluctuations in iron-based compounds |9|) yields 
Tc = 5 Kelvin consistent with the experimental transi¬ 
tion temperature of 8 Kelvin in FeSe [3l |4] (note that 
allowing g{k) to be a sum of symmetry functions further 
increases A m and our calculation therefore provides a 
lower bound for the transition temperature). In the cur¬ 
rent work, we have not considered the effect of phonons 
which was recently shown to be significant in the FeSe 
monolayer [39| . Future work is needed to investigate the 
interplay of spin fluctuations and phonons in iron-based 
superconductors. 
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Note added: After submission of the present 
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mediated superconductivity in iron selenide was pub¬ 
lished [40] . 
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